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Abstract 

The effects of isospin violation on the neutral weak magnetic form factor of the proton 
are studied using two-flavour chiral perturbation theory. The first nonzero contributions 
appear at 0(p 4 ) in the small-momentum expansion, and the 0(p 5 ) corrections are also 
calculated. The leading contributions from an explicit A (1232) isomultiplet are included 
as well. At such a high order in the chiral expansion, one might have expected a large 
number of unknown parameters to contribute. However, it is found that no unknown 
parameters can appear within loop diagrams, and a single tree-level counterterm at 0(p A ) 
is sufficient to absorb all divergences. The momentum dependence of the neutral weak 
magnetic form factor is not affected by this counterterm. 



1 INTRODUCTION 



The first measurement of the proton's neutral weak magnetic form factor has been re- 
ported recently by the SAMPLE Collaboration at MIT/BatesQ; subsequently, a series of 
precision experiments has gotten underway at the Jefferson Lab [[J. Assuming only two 
quark flavours with exact isospin symmetry and neglecting electroweak radiative correc- 
tions, the weak form factor can be expressed in terms of the familiar electromagnetic form 
factors as follows: 

G p m Z (q 2 ) = \{G\M 2 ) ~ GIM 2 )] - G\ 4 {q 2 W0w ■ (1) 

Electroweak radiative corrections have been discussed in Ref. ||. 

There is presently a great deal of interest in determining the contribution due to 
strange quarks which simply appears as a new term added to the right-hand side of 
Eq. flU). It is important to notice that isospin- violating effects also appear as a new term on 
the right-hand side, even in the absence of strange quarks. There have been some attempts 
to estimate the isospin- violating effects by using constituent quark models f| §, [7| and a 
light-cone meson-baryon fluctuation model[f|. 

In the present work, the effects of isospin violation are studied using two-flavour heavy 
baryon chiral perturbation theory (HBChPT) |9|, [L(J], which is nicely suited to the task. At 
small momentum transfer, HBChPT is a systematic expansion in small parameters with 
dynamics arising from the propagation of pions and photons in the presence of a single 
baryon. The remaining short-distance physics can only appear as low-energy constants 
(i.e. parameters in the HBChPT Lagrangian). The spontaneously-broken chiral symmetry 
of QCD is respected in HBChPT by construction, and the explicit breaking due to current 
quark masses can also be included in a systematic way. 

Isospin violation occurs in nature through electromagnetic as well as strong interac- 
tions, so it is necessary to include the effects of virtual photons in HBChPT. The required 
Lagrangian has been constructed by Meissner and Steininger||l 1||. 

In principle, t/iree-flavour HBChPT can be used to explicitly include strangeness in 
the meson-cloud contribution to the neutral weak form factors, and some work in this 



direction has recently been reported. [12| However, the chiral expansion is not as well- 
behaved for strange quarks as it is for up and down quarks, and since the present work 
is only concerned with isospin-violating effects it is preferable to work with two-flavour 
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HBChPT. 

The neutral weak vector form factors of the nucleons are defined here in the notation 
of Dmitrasinovic and Pollock||, 

(N(p + <f)\±(uj»u-d^d)\N(^ = u{p + q)^-r d Ff +n ± u - d Fr n )l, 

u(p), (2) 



,}_/u~d rpp+n I u-d Tpp-n\ ^JlvQ 

2 l 2 2 1 2M N 



N(p + q)\^(u-f^u + d-f^d)\N(p)) = //(/7+ <i) 



_/u+dj?p+n _j_ u+d pp-n\ 



,}_/u+dpP+n I u+dpp-n\ WjlvQ 

2 l 2 2 ' 2M N 

with M^v denoting a nucleon mass. Walecka-Sachs form factors are defined by 



nip), (3) 



l GW) = l F[(q 2 ) + l Fi(q 2 ), (5) 



l GW) = + (4) 



where i = u ± d and j — p ± n. When isospin conservation is not enforced, Eq. ([!]) 
generalizes to 

(g 2 ) = \[G\M 2 ) - GlM 2 )} - G\M 2 We w - ^(g 2 ), (6) 

where 

G M (g2) s -Wc^V) _ (7) 

is the isospin-violating term. 

In section 0, it is shown that isospin-violating contributions do not appear up to and 
including 0(p 3 ) in HBChPT. Section [3] presents and discusses the 0(p A ) calculation. At 
this order, the isospin-violating % F{ form factors contain no unknown parameters. Each 
% Fl requires a single counterterm, but the momentum dependence is not affected by the 
counterterm. 

In section |J the calculation of G u ^{q 2 ) is extended to next-to-leading order, 0(p 5 ). 
In principle, a large number of low-energy constants could appear, some within loop 
diagrams and others as tree-level counterterms, but it is shown that the loop integrals 
are finite and no new low-energy constants appear at this order. The ratio of next-to- 
leading versus leading order contributions provides some indication of the behaviour of 
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the HBChPT expansion. For the derivative of G^f(q 2 ) at q 2 = 0, this ratio is close to 
1/2. 

Section [5] evaluates and discusses the contributions to G^f(q 2 ) that arise when the 
A(1232) is included explicitly in the chiral Lagrangian. Section [| offers a determination 
of the pion-cloud contribution to the sole remaining parameter in G^f(q 2 ), and then 
summarizes the complete HBChPT result for G^f(q 2 ). 

2 NO CONTRIBUTION UP TO 0{p 3 ) 

The Lagrangian of HBChPT is written in the form 

r - /•(*) i r( 2 ) . r( 3 ) . r( 4 ) . r( 5 ) . (*\ 

where the superscripts denote powers in the "momentum" expansion, which is actually 
a combined expansion in various small quantities. A covariant derivative counts as one 
power (therefore the field strength of an external current counts as two powers) and a 
current quark mass counts as two powers (recall that m 2 ~ m q ) . These dimensionful quan- 
tities are small relative to both the chiral scale and the nucleon masses, AirF^ ~ M p M n . 
Eq. (Rp represents an expansion in the inverse nucleon masses as well as the chiral ex- 
pansion. Furthermore, the effects of virtual photons can be organized according to an 
expansion in the electromagnetic coupling, and it is convenient to use 0(e) ~ Q(p)||ll]|, 
which allows the virtual photon effects to also be incorporated into the generalized mo- 
mentum expansion of Eq. In the present work, u O(p n )" is used to denote n powers 
of any of these small quantities. 
The lowest-order Lagrangian is 



tf$f = N v (ivV + g A S-u)N v , (9) 



where 



N v (x) = exp[iM vx]hl + f)il)(x), (10) 

S H = 7fls0^ u V V , (11) 



Ufl = iu ] (d^-ir^)u-iu(d^-U^)u\ (12) 
V M = 0„ + r„-i«W, (13) 
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(14) 



and Mq is the lowest-order nucleoli mass. External vector and axial vector fields are 



included via r. 



+ A^ and 



— Ap, and u is a nonlinear representation of the 



pion fields, for example 



u 



exp 




71 



-It 



(15) 



The parameter F corresponds to the pion decay constant in the chiral limit (normalized 
according to F n w 93 MeV). 

The following relations will also be useful: 



S-v 



0. 

bKz pVpuJ u u ■ 



(16) 

(17) 
(18) 



The relativistic currents required for this work can be re-expressed as a 1/Mn expansion 
between HBChPT spinors. In the rest frame of the initial nucleon, one finds 



2M 



N„ 



V » + 2M 



1 



N 



M, 



-le 



N 



pvpa 



<?v p S a + O 



1 



M 2 N , 



N 



-ie. 



q u v p S w + 



+ 



(19) 
(20) 



where = (1, 0, 0, 0), is the 4-momentum of the incoming vector current and Mjy is 
the physical mass of the nucleon. 

It is a simple matter to determine from C^n the tree-level coupling of an external 
neutral vector field to a nucleon. The result is 



u—dpp—n u+dpp+n 

u—dpp+n u+dpp-n 



l Fi 



0, for all 



(21) 
(22) 



No contributions to the % F^ form factors can appear at this order due to the explicit factor 
of q/M^ in Eqs. @ and (§). Also, loop graphs are forbidden at 0(p) by the standard 
power counting of HBChPT, which assures that one-loop diagrams constructed from C^ N 
begin at 0(p 3 ). 
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The isospin-violating form factors must contain m u — or a virtual photon (with 
the associated factor of e 2 ), and are therefore suppressed by at least two powers relative 
to the isospin-conserving ones. This guarantees that u ~ d pv +n an d u + d F^~ n must remain 
zero until 0{p i ). One might expect u - d Ff +n and u+d Ff~ n to be nonzero at 0(p 3 ), but 
in fact they also remain zero until 0(p 4 ) due to Noether's theorem and the fact that the 
vector currents uj^u and d'y^d are each conserved in QCD and QED (recall that weak 
radiative corrections are being neglected). Noether's theorem requires that the l F{ form 
factors of Eqs. (j2lD and ([22|) do not get renormalized at q 2 = 0, but explicit factors of q 2 
cannot appear until 0(p A ) in loop diagrams or 0(p 5 ) at tree level. 

Finally then, it is concluded that all isospin-violating form factors remain zero up to 
and including 0(p 3 ). It is instructive to verify this fact by a direct calculation using the 
chiral Lagrangian which has been written down in its entirety at this order. [|I3|, TA, Tl\ 



Following the notation of Ref. |13|] where field redefinitions have been used to remove 
"equation of motion" terms from the Lagrangian, the 0(p 2 ) terms which affect the neutral 
vector form factors are 



xri 2 ) — AT 
0J -ixN — ly v 



where 



N v , (23) 



Unv = u{d/ u - dj^-i^^4 u })v) + v) {d^r v - d v r^-i[r^r v ])u, (24) 
«J2 = dpvP-drvU. (25) 

Adding these terms to the lowest-order Lagrangian of Eq. @ leads to the following 
coupling of an isoscalar vector current to a proton: 



N v . (26) 



According to Eq. (0), the expression for an isovector current is obtained by the replace- 
ment a-j 2a 6 . q^ is the incoming 4- momentum of the vector current, and the 4- vectors 
and hp are defined by 

p„ = MoVp + (27) 
where p^ is the 4-momentum of the incoming nucleon. 
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One is free to work in the rest frame of the initial nucleon, p^ = (Mat, 0, 0, 0), and to 
choose t> M = (1, 0, 0, 0). For on-shell nucleons, these choices imply 



= (M N -M ) VfM , (2i 
2M N 



V -Q = (29) 



and lead to 



u+drpP-n u-drpP+n (M n M p ) 

Fl = Fl = Mo ' (30) 

u+dpP-n = u-dpp+n = Q ^ 



Recalling that the neutron-proton mass difference is 0(1/M O )[|I3], [TJJ, one concludes that 
the nonzero result of Eq. (^) is suppressed by two powers of 1/M relative to the leading 
isosp in-conserving result. In other words Eq. (|30D contributes at 0(p 3 ), and therefore the 
isospin- violating form factors all vanish up to 0(p 2 ), as expected. 

When the calculation is extended to 0(p 3 ), Eq. (|30|) is not the only contribution. C^, 
which includes virtual photons as well as strong interactions, contains approximately 40 
new parameters. |T3], |T4] , |TT| However, none of these parameters can contribute to the 



isospin- violating form factors at 0(p ), and the only terms that do contribute are 



6Kn ~ 2hlf v 



1\ . / 1 
4 



fl 6 - Z ) + ( a 7 - 7 ) V »v 



e^ pa S a iV p + h.c. N v . (32) 



A simple calculation shows that these terms do not affect l F[, and their contributions to 
% Fl are effectively 0(p A ), so they are negligible in an 0(p 3 ) calculation. 

In a general HBChPT calculation, one-loop diagrams built from interactions can 
contribute at 0(p 3 ), but none exist which can contribute to the isospin- violating form 
factors of interest here. However, there is a contribution from wave function renormaliz- 
ation[l5|, 



and it precisely cancels the 0(p 3 ) effect found in Eq. ([II]). All 0(p 3 ) effects have now 
been discussed, so the isospin-violating form factors do indeed remain zero at this order. 
It is easy to verify that the final results of this section remain unchanged if one does 



not employ the field redefinitions of Ref . |I3[ , even though unphysical intermediate steps 
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may differ. For example, the wave function renormalization constant at 0(p 3 ) becomes 
independent of the nucleon mass when the field redefinitions are not used, so Z p = Z n , 
but this change is exactly compensated by the extra equation-of-motion term that would 
be present in £^yfl3||, 

^N^^i + ^fNvivVfNv. (34) 



3 LEADING ORDER, 0(p 4 ) 

Having verified explicitly that the isospin-violating neutral weak vector form factors are 
exactly zero up to 0(p 3 ), the calculation will now be extended to 0(p A ) where a nonzero 
result does exist. The complete Lagrangian £^y has not been written in the literature (but 
an effort is underway: the counterterms required for renormalization have recently been 
listed 0]). However, the contribution of C^ N to the present study is a simple constant and 
the full Lagrangian need not be constructed here. Of greater interest are the g 2 -dependent 
0(p 4 ) effects which come from loop diagrams constructed using + C^N- 

All of the pion-loop diagrams which contribute to the isospin-violating form factors at 
0(p A ) are shown in Fig. |IJ (Diagrams with virtual photon loops will be discussed later 
in this section.) Each of them will be evaluated in the rest frame of the initial nucleon 
with the "velocity" parameter fixed to = (1, 0, 0, 0). Final results for the form factors 
do not depend on these choices. 

The pion propagator and the vector-pion couplings are obtained from the lowest-order 
chiral Lagrangian for mesons, 



£l 2 2 = -Tr 



D^D^U + 2B 





where U = u 2 and 



4 

+e 2 C I " I / I J P . (35) 

^ -1/3 ) \ -1/3 j V ; 

D tl U = d fl U-tr tl U + tU£ fl . (36) 



Notice that the physical pion masses are nonzero for two reasons: current quark mass 
effects (the parameter B in C^l) and electromagnetic effects (the parameter C in C^). 
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Figure 1: Contributions to the isospin- violating vector form factors of a nucleon at 0(p 4 ). 
A dashed line represents the sum over charged and neutral pions. The solid dot denotes an 
insertion from £^ N + £^y, a cross denotes an insertion from C^ N , and all other Feynman rules 
come from C^h + C^l ■ 



The wave function renormalization constant for C % n with the field redefinitions of 
Ref. [13[] has been determined previously up to 0(p 3 )|T5], [T(J, but the present work will 
require an extension of the isospin- violating part to 0(p 4 ). The relevant diagrams are 
displayed in Fig. |[ Using dimensional regularization, the calculation of a diagram without 
insertion proceeds as follows: 

d d £ ( i N 



an 0{p 2 
a A z 



A-d 



(2vr) d \£ 2 - mK + ie 



4(4ttF)< 



v-k ( mlo - ~(v-k) 2 



9a 
' F 

2 



S- 



+~v ■ /,' {mlo 



[vkf 



4/ 2 



d 

{vkf 



1 



v ■ (k -t) +ie) 
- 7 + ln(47r) - In 



9 a 



S ■£) (37) 



3/2 / 7T V-k \ 

— h arcsm 

\ 2 m^o / 



(38) 



Although this expression contains no explicit quark masses, it can still violate isospin by 
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A z ~ -l — ' A z - ' X ' 

+ + 

A z"^ — ' A z ~ ' X 

Figure 2: Contributions to the isospin- violating piece of the nucleon's wave function renormal- 

(2) 

ization constant. A cross denotes an insertion from and all other Feynman rules come 

from C^jf + CnJ ■ 

virtue of the on-shell relation between v-k and the nucleon masses as given in Eq. fl2*8p. 
The wave function renormalization constant is defined to be the residue of the nucleon 
propagator at the on-shell point (multiplied by i), so Eq. ( |3~8D should be viewed as a 
function of x = v-k — + M . The resulting contribution of this one diagram to the 
difference between proton and neutron wave function renormalization constants is 

5(Z p - Z n ) = i± [ a A z (p) - a A z (n)} x=Q = -|(M B - M p )j^ + 0(M n - M p f. (39) 

Similar contributions are made by the other diagrams in Fig. E3, 



. d 

dx 


~ b A z ( P ) 


- b A z (n)_ 


x=0 


~(M n - 


. d 

dx 




- C Az(n)] x=0 = 


-3(M n 


. d 

dx 


' d Az( P ) 


- d A z (n) 


x=0 


-3(M„ 



+0(M n -M p f, (40) 



(4ttF) 



,2 , 



(ArrFf 



+ 0(M n -M p y, (41) 



ir&m^ +0{Mn _ Mp) 2_ (42) 



PJ (4ttF) 2 

Adding these four contributions to the lower-order result of Eq. (1331) yields 



^ - z " = - 6(M » - M ^74^r • (43) 

The expressions for the matrix elements corresponding to the diagrams of Fig. [I] are 
now presented in the rest frame of the initial nucleon with = (1,0,0,0), on-shell 
external nucleons, and an isoscalar vector current with incoming momentum 

a A^ +d \p)- a A^ +d \n) = const x i(M n - M p )it^ pa q u v p S% (44) 
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h A^ +d \p)- h A^ +d \n) 

c A^ +d \p)- c A^ +d \n) 
l A^ +d \p)- d A^ +d \n) 



-i(M n - M p 



6v 



M (4ttF) 2 
(2m 7r H 



k9a 



f A^ +d \p) - f A^+ d \n) 
Their sum is 

A^ +d \p) - A^ +d \n) = i(Af, 



p; (4ttF) 2 ' 
A( u+d )(p)- e 4«+ d )(n) 

2(2771,,+ -77V)), 

9 A^ +d \p)- 9 A^ +d \n) = 0. 



-u; M (M n -M„)-^- 
4 MV n W (4vrF)< 



2M 2 



+ const x ien Upa q v v p S a 



(46) 

(47) 
(48) 

(49) 



The contribution from Eq. fl32|) has also been absorbed into the unspecified constant. 
According to Eqs. ([19]) and (p0|), the term containing ie flupa q u v p S a in Eq. (|49| ) is the 
leading contribution to " +d F2' _n (g 2 ). If there had been a term containing v^, it would 
have been the leading contribution to u+d F p ~ n (q 2 ). As discussed in the previous section, 
u+d F p ~ n (0) = is required by Noether's theorem, and this is explicitly verified by the 
exact cancellation of terms proportional to in Eqs. (44)-(f48|). 

There is also a term in Eq. ([fjj) that is proportional to q p , and its origin is more 
subtle. Instead of contributing to the isospin violation of the form factors, this term 
is a consequence of the nucleon mass dependence in the currents themselves. Recall that 
Eqs. (|i~9|) and ( p0| ) are expressed in terms of the physical nucleon mass rather than the 
bare mass. Thus, 



(M n - M 3 



2M 2 



q» + 2ie ilvpa q v v p S a + O 



M \ 



N v . 
(50) 

while the term containing 



This term containing is precisely the one shown in Eq. (| 
the antisymmetric tensor is absorbed into the unspecified constant in Eq. ([49]) . Therefore, 
the 0(p 4 ) expressions for the "77 + d" vector form factors are 

u+d Ff- n (q 2 ) = 0, (51) 



u+d T-ip—n, 
^2 



\q 2 ) = const = u+d K p ~ n . (52) 

The form factors for an isovector ( "77 — d" ) vector current are less trivial because they 
receive a nonvanishing contribution from those diagrams in Fig. [I] where the vector current 
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couples directly to a pion. Those are the diagrams which produce momentum dependence 
at 0(p 4 ). The complete set of u u — d" contributions leads to 



^Ff+V) = -T^rn« + {M n -M p) 




\ 



1 — X(l — X 



2 

m + 



- d Ft n (q 2 ) = u ~ d ^ +n + W f A ^ (M n - M p ) f 1 dx In ( 1 - g(l - g)-^-) 



(53) 



(54) 



The numerical value of u - d ^p +n is not specified by chiral symmetry. It should be noted 
that Eq. ([53]) satisfies the constraint that "~ d Ff +n (0) = 0, which is a nontrivial check of 
the algebra since it results from a cancellation among the diagrams of Fig. [I]. 

Throughout this section, pion-loop diagrams have been evaluated while the effects 
of virtual photon loops were tacitly omitted. This omission can be justified by simple 
power-counting arguments. To begin, recall that the chiral Lagrangian is order-by-order 
renormalizable, and that any two-loop diagram cannot contribute before 0(p 5 ). (Photon 
loops have the same power counting as pion loops due to the association 0(e) ~ 0(p) JITJ.) 
Up to 0(p A ) then, the only effects of virtual photons come from the addition of a single 
virtual photon to the simplest tree-level form factor diagram. 

Such a one-photon- loop diagram could offer an 0(p 3 ) contribution to l F{(q 2 ). However, 
this form factor is required to vanish at q 2 = 0, so there must be an extra factor of q 2 /Mq at 
least, which then contributes at 0(p b ). This stands in contrast to the pion-loop diagrams 
which contain the extra mass scale m n ] the ratio q 2 /m 2 does not lead to extra suppression 
in the HBChPT expansion. 

A one-photon-loop diagram could also offer a contribution to *iV(g 2 ), which begins 
at 0(p 4 ) rather than 0(p 3 ), due to the explicit factor of q/M^ in the definition of these 
form factors, Eqs. (fj) and (Q). The contribution must be a simple constant because any 
q 2 dependence would require extra factors of 1/M and would thus contribute at a higher 
order. The constant 0(p A ) contribution from virtual photons simply adjusts the numerical 
value of the unspecified parameters, % k? in Eqs. (|52j) and (|54|). 

The 0{p A ) isospin- violating u u — d n form factors are plotted in Fig. § for < — q 2 < 
0.25 GeV 2 , with «-<V +n set to zero and the numerical values of all other quantities set 
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.030 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




.00 .05 .10 .15 .20 .25 

-q 2 [GeV 2 ] 

Figure 3: Leading-order HBChPT results for the isospin-violating isovector form factors. For 
this plot, the unspecified constant in u ~ d pP +n has been chosen such that the form factor vanishes 
at q 2 = 0. 

to those of the Particle Data Group. [|18| Uncertainties are not shown in Fig. [| As will be 
discussed in the following section, the dominant uncertainty comes from the truncation 
of the HBChPT expansion at 0(p 4 ), and an estimate of this uncertainty will be obtained 
from the 0(p 5 ) calculation. 

Independent of the value of u ~ d Ki P+ n j Fig. |3] indicates that both ll u — d" form factors 
are monotonically increasing at leading order, with u ~ d pP +n increasing more quickly than 

According to Eqs. it is the sum of these two form factors which is relevant 

to the proton's neutral weak magnetic form factor. However, l F( and differ by an 
explicit power of Mn due to their definition, Eqs. (||) and (||). Thus, the leading order 
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result from HBChPT is 

G u A :f(q 2 ) = - - l 4f^{M n - M p ) I' do; In (l - x(l - x)-^-) . (55) 

V*xr y Jo \ m w+ J 

The 0(p 4 ) result for «~ rf Ff +n is a next-to-leading-order correction to Gfy , and it is reas- 
suring to see from Fig. |3]that its g 2 -dependence is smaller than the leading g 2 -dependence. 

The remaining next-to-leading-order effects come from an 0(p 5 ) calculation of l F%, 
and are discussed in the following section. 

4 NEXT-TO-LEADING ORDER, 0{p 5 ) 

The next-to-leading corrections for both the % F{ and l F% isospin-violating form factors 
occur at 0(p 5 ) in HBChPT. The goal of the present section is to complete the next-to- 
leading calculation of G^f. l F( is itself a subleading contribution to G^f, so the 0(p 4 ) 
result of the previous section is a next-to-leading order effect. Therefore only needs 
to be calculated at 0(p 5 ). 

To construct the set of contributing Feynman diagrams, consider first those diagrams 
which contain no propagating pions (i.e. tree-level diagrams or diagrams containing any 
number of photon loops). The only available dimensionful parameters are the nucleon 
mass and momentum transfer, and it is easily concluded that such diagrams cannot con- 
tribute to the isospin-violating % Fl at 0(p 5 ). Also, no contribution emerges from any 
diagram which contains both one photon loop and one pion loop. 

A diagram containing two pion loops can only contribute at 0(p 5 ) if all Feynman 
rules come from C^ N + C^. However, no isospin violation is contained within C^ni an d 
the m 7T + — m^o mass difference does not affect the form factors of interest here, so only 
one-pion-loop diagrams can contribute. 

The set of diagrams which do contribute includes b A fl , ? and 9 in Fig. [I], plus the 
diagrams of Fig. ^. Although contains about 40 parameters that are unconstrained 
by chiral symmetry and could in principle appear within 0(p 5 ) loop diagrams, none of 
them contribute to this calculation. Dimensional arguments do not permit a tree-level 
0(p 5 ) counterterm for % Fl, (essentially because the small expansion parameters without 
uncontracted Lorentz indices tend to come in pairs at tree-level, such as m 2 , e 2 , or q 2 , 
but 0(p 5 ) would require an odd power), so the total loop calculation must be finite. 
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Figure 4: Contributions to the isospin-violating form factors that begin at 0(p 5 ). A dashed 

(2) 

line represents the sum over charged and neutral pions. A cross denotes an insertion from C ^, 
and all other Feynman rules come from 4at + ^1 



Summing the 0(p 5 ) contributions and adding them to the 0(p A ) results from the 
preceding section produces the full isospin-violating contribution to the proton's neutral 
weak magnetic moment, as computed within HBChPT (without explicit A(1232) fields), 



d/ 2\ 



u+d ^p-n u—d ^p+n 



i6g\M N 

(4ttF) 2 



(M n - M v ) C dxln ( l-xil-x)-^—] 

Jo \ 171*+ J 



[ n p) (4ttF)2 



+ 



1 r 1 



6 Jo 



dx I 1 — x(l — x) 



/i N 2 



xd — x) 



m 



(56) 



The parameters ciq and 07 have been re-expressed as the nucleon magnetic moments via 
Eq. 



«6 
7 



4/Uat 

/^p "I - /^n 
2jl N 



1.176, 
0.440, 



(57) 
(58) 



plus higher-order corrections that are not needed for the present work, /i^r is the nuclear 
magneton. 
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Figure 5: The leading order (LO) and next-to-leading order (NLO) HBChPT results for the 
isospin-violating contribution to the proton's neutral weak magnetic moment. For this plot, the 
unspecified constant at LO has been chosen such that the form factor vanishes at q 2 = 0. There 
are no unspecified parameters in the NLO contribution. 

Fig. [5] contains a plot of G^f(q 2 ) versus q 2 at leading order (LO) and next-to-leading 
order (NLO), with the only unspecified quantity u+d ^p^ n — u~d K p+n ge ^. ^ Q zero ^he NLO 
effects serve to soften the g 2 -dependence of the form factor. Fig. [5] indicates that the 
NLO corrections to G^f(q 2 ) total roughly 0.01 at q 2 = 0, and grow to about 0.02 near 
q 2 = 0.1 GeV 2 . 

To determine how well the HBChPT expansion is working for this observable, it is 
useful to consider the derivative of the form factor at q 2 = 0, 



d(-q 2 ) 



G U M d (0) 



8g 2 A M N (M n -M p ) 
3m 2 + (4vrF) 2 



97rm 7I 



8M 



N 



(59) 



Thus the ratio of magnitudes of the NLO/LO contributions is 9nm n + /8Mjy ~ 1/2. Taking 
this as representative, the uncertainty from the neglect of NNLO effects could be roughly 
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half of the NLO contribution. These uncertainties will be discussed further in section |6]. 



INCLUDING THE DELTA RESONANCE 



Conspicuous by its absence from the preceding discussion is the A(1232) isomultiplet. 
The effects of an infinitely-heavy A(1232) are accommodated within the numerical values 
of the HBChPT parameters, but the rather small A-iV mass difference observed in nature 
raises the possibility of substantial corrections to the Ma oo limit. 

The incorporation of an explicit A(1232) field into HBChPT was initiated by Jenkins 
and ManoharH and has been employed by various authors. [19 , |20|j In this section, the 



formalism developed by Hemmert, Holstein and Kambor in Ref. |2(J is used to calculate 
the leading A(1232) contributions to the isosp in- violating l Fl and form factors. 

As is familiar from relativistic approaches to spin-3/2 field theory, the Lagrangian 
contains a vector-spinor and one then employs a projection operator to isolate the spin- 
3/2 piece. For example, the lowest-order propagator in "<i" spacetime dimensions is 

4 



v -k — A + ie 



d- 1 



(60) 



where the incoming 4-momentum is M& v^ + Ma,o is the lowest-order A(1232) mass, 
and 

A = M A ,o - M . (61) 

As expected, the leading contributions of the A(1232) to the isospin-violating form 
factors appear at 0(p A ). The relevant diagrams are displayed in Fig. [?], and the terms 
required from the leading-order Lagrangian are 



<5£ (1) 



7rA 



nNA 



_rpH 

i 

dirNA 



iv-D lJ - 5 lJ A 



T 



A* J' 

T?u]N + NufiT? 



where is the vector-spinor, is defined by 

< = ^Tr (r^) > 
(t 1 is a Pauli matrix in isospin space) and the covariant derivative is 



r,-4 s) ) 



2ie»J"r 



fc-pfc 



(62) 
(63) 

(64) 
(65) 
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Figure 6: Contributions of the A(1232) to the isospin-violating vector form factors of a nucleon 
at 0(p 4 ). A dashed line represents the sum over charged and neutral pions, and a double line 

(2) 

represents the A(1232). A cross denotes an insertion from CjlL, and all other Feynman rules 



come from + C^ NA 



+ £ (2) 



The only insertions from that contribute are the A(1232) mass corrections, which 
arise from electromagnetic as well as strong interaction effects. 

The authors of Ref. |2(| chose not to perform the A(1232) field transformation that 
would have removed "equation of motion" terms from C^. None of the calculations in 
this section depend upon whether or not the transformation is performed. 

The leading contributions of the A(1232) to nucleon wave function renormalization are 
shown diagrammatically in Fig. 0. When their contribution is added to the non-A(1232) 
result of Eq. fl43|) , the full isospin- violation due to wave function renormalization is found 
to be 



Mn - Mp -6(M n - M y 9 ^: 

M \ n p) ( 47rF ) 2 
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Figure 7: Contributions of the A(1232) to the isospin-violating piece of the nucleon's wave 
function renormalization constant. A double line represents the A(1232). A cross denotes an 
insertion from C%\, and all other Feynman rules come from + C^NA "F 



where 
R{m 2 ,q 2 ) 



--[4(M n - M p ) + (M A o - M A+ ) - 3(M A - - M A++ )]R(ml + , 0) 
16 

-—[M n -M p - (M A o - M A +))R(m 2 , 0). 



(66) 



A 2 

(4ttF) 2 

a r 1 



A f - - 7 + ln(47r) - ln^ ] - A f dx In ( 1 - x(l - x)-^] 

\e fi z J Jo \ m 2 J 

x ( x{l-x)q 2 \ 
\A 2 — m 2 + x(l — x)q 2 J 

j 1 I 2A 2 — m 2 + x(l — x)q 2 x(l — x)q 2 [m 2 — x(l — x)q 2 } 

A 2 - m 2 + x(l - x)q 2 3[A 2 -m 2 + x(l -x)g 2 ] 3 / 2 

A 



x In 



^Jm 2 - x{\- x)q 2 \m 2 - x(l~ x)q 



A 2 



- 1 



(67) 



The set of 0(p A ) diagrams, contained in Fig. || produces the following contributions 
to the isospin-violating form factors: 



u+dpp-n/2 



s 



u+d F%~ n {q 2 







(68) 
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- d F[ +n (q 2 



•y [2(M n 



(M A o - M A+ ) - 3(M A - - M A ++)] 



x[i?(m 2 + ,g 2 )-#(m 2 + ,i 



r Mjv [2(M n - M p ) - (M A o - M A 



(69) 

3(M A - -M A++ )] 




(70) 



Notice that the total contribution made by the A(1232) loop graphs to each of the "w + <i" 
form factors exactly vanishes. For the isovector ( u u—d" ) current, the various contributions 



to 



u-d pP+n 



F{ add in such a way that all divergences cancel and the familiar constraint, 
u ~ d F[ +n (0) = is satisfied. As expected, the contribution to u - d FP +n [ s not finite, and 
the full Lagrangian contains an 0(p 4 ) counterterm which absorbs this divergence. 

Plotted in Fig. |§| are the leading results, i.e. the 0(p 4 ) results, for the two "w — d" 
form factors with and without the explicit A(1232) contribution, and with the unspecified 
constant subtracted from u ~ d F$ +n . The effects of the A(1232) are significantly smaller 
than the NLO corrections of Fig. [5]. The value g n NA = 1-05, recommended in Ref. [ B0| , 
has been used along with the following mass differences: 



A = 0.293 GeV, 
M A o - M A+ = -(M A - - M A++ ) = .0013 GeV. 



(71) 
(72) 



There are sizable experimental uncertainties on these inputs, but no reasonable choices 



can make the A(1232) contribution to G^f(q 2 ) grow larger than the NLO correction of 
Fig. |. 



6 DISCUSSION 

The proton's neutral weak magnetic form factor G p ^f is of great interest, both exper- 
imentally and theoretically, because it is sensitive to strangeness within the nucleon. 
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Figure 8: Leading-order HBChPT results for the isospin-violating isovector form factors with- 
out (solid lines) and with (dashed lines) the explicit A(1232) isomultiplet. For this plot, the 
unspecified constant in u ~ d F^ +n has been chosen such that the form factor vanishes at q 2 = 0. 

Experiments actually measure the sum of strangeness and isospin-violating contributions 
(labeled G S M and G^f respectively), so it is advantageous to understand the effects of 
isospin violation as thoroughly as possible. 

In this work, heavy baryon chiral perturbation theory (HBChPT) has been used to 
study isospin violation in the absence of strange quarks. HBChPT is an expansion in 
momenta that are small compared to the chiral scale, A x ~ m p ~ Mjy ~ Virtual 
photons can also been included according to the usual q;qed expansion. 

First, G^f(q 2 ) was calculated at leading order (LO) and next-to-leading order (NLO) 
for a single nucleon surrounded by a cloud of pions and photons. The photon cloud 
was found to contain no momentum dependence at this order, in contrast to the pion 
cloud. The ratio of NLO/LO contributions was used as an estimator of the systematic 
uncertainty of the HBChPT expansion. Then the contribution of the A(1232) isomultiplet 
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was evaluated and found to be negligible in comparison with the systematic uncertainty. 

The final expression is given by Eq. (]56f) in terms of a single unspecified parameter, 
u+d K p-n _ u-d^p+n^ ^his parameter contains physics of two types: a low-energy contri- 
bution from the pion cloud, and a higher-energy contribution which is unconstrained by 
HBChPT. 

The pion-cloud contribution to the unspecified parameter is easily determined by re- 
doing the HBChPT calculation with a momentum cutoff, A, instead of dimensional reg- 
ularization. This cutoff represents the separation scale between the "low-energy" and 
"higher-energy" regions. In principle, the choice of A does not affect G^f(q 2 ) since only 
the sum of low-energy and higher-energy pieces is relevant, but in HBChPT the higher- 
energy piece is undetermined so it is preferable to make a physical choice for A. 

Clearly A cannot be larger than A x , since HBChPT fails above this scale. In a series 



of recent papers|21]], Donoghue, Holstein and Borasoy have argued for an HBChPT cutoff 
that is not too far above 1 fm _1 « 200 MeV, corresponding to the measured size of a 
baryon. Above this approximate scale the substructure of a nucleon can begin to be 
relevant, but is not accurately represented in HBChPT. 

The only divergence in the non-A(1232) piece of G u £^{q 2 ) comes from a single integral 
which appears in diagrams f and 9 A fl of Fi gure |TJ. When dimensional regularization 
is replaced by a simple cutoff for the momentum integral, the following relationship is 
obtained: 

' ' -T + lnM-ln^l - fln4__?V (73) 



4 — d /i 2 I \ m 



-+ 



Choosing A = 400 MeV, the low-energy contribution of the pion cloud is found to be 

0.014. (74) 



u+d p—n u-d J)+ij 



pion cloud 



Notice that the dependence on A is only logarithmic, so the result is not overly sensitive 
to the chosen numerical value of the cutoff. The A(1232) contribution also contains a 
logarithmic dependence on A, but the resulting pion-cloud contribution is negligible in 
comparison to the uncertainties coming from the HBChPT expansion (recall section ^). 

The remaining contribution to the G u ^f counterterm is the "higher-energy" contribu- 
tion. It is unspecified in HBChPT by definition, and a precise numerical prediction is 
therefore beyond the scope of this work. If the HBChPT expansion is to be well-behaved, 
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then the higher-energy contributions must respect the established power counting. Per- 
haps the most obvious examples of higher-energy physics are the p and u vector mesons. 
Simple power counting estimates for the tree-level vector meson dominance diagrams 
indicate that their leading contribution resembles 

M N (M n -M p ) 



OC 



p,u> 



m 



(75) 



p,LU 



which is 0(p A ), as required by HBChPT. 
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Figure 9: The HBChPT prediction for the isospin-violating pion-cloud contribution to the 
proton's neutral weak magnetic form factor. The two solid lines represent the central value with 
and without a subtraction at q 2 = 0, and a pair of uncertainty estimates are provided for each 
solid line according to Eq. (|Tq). The logarithmic divergence is cut off at 400 MeV, as discussed 
in the text. 

The main conclusion of this work is summarized by Fig. ^ which shows the full 
contribution of the pion cloud to (g 2 ), for < -q 2 < 0.25 GeV 2 , up to next-to- 
leading order in the HBChPT expansion. To aid a discussion of g 2 -dependence, the same 
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quantity is shown after G^f(0) has been subtracted. In each case, a pair of uncertainty 
bands is shown, representing two estimates of the error associated with the neglect of 
NNLO contributions, 

I i|NLO| , wideband , , 

|NNLO| ~ < 2 1 , , . (76) 

I -^-|NLO| , narrow band 

The wide-band error estimate purports that the ratio of 1/2, taken from |NLO|/|LO| in 
Eq. (|59|) for the derivative of G^f, might be a reasonable indicator of the NNLO uncer- 
tainty. The narrow-band error estimate employs a generic HBChPT expansion parameter. 

Fig. |9] indicates that the g 2 -dependence, as determined by G^f(q 2 ) — G^f(0), occurs 
on the scale of a few times 0.001, but is typically less than 0.01. This result is easily 
understood: HBChPT demands that the momentum dependence of this isospin-violating 
form factor is proportional to (M n — M p ) / ~ 0.001. Recall that no contributions from 
171^+ — m^o or any of the HBChPT parameters were found. 

However, Fig. |3| also indicates that the contributions of the pion cloud to G^f(0) are on 
the scale of a few times 0.01. The origin of isospin violation is still solely (M n — M p )/Mn ~ 
0.001, but the numerical coefficients are larger than those for the g 2 -dependence. 

It is interesting to compare the results of Fig. |] to the findings of other authors. 
Dmitrasinovic and Pollock have used a nonrelativistic constituent quark model to find [[J 

" +d GV(0)l = \ u - d GT(0)] - 0.008 =► [Gtf(0)l = 0. (77) 



DmiP 



DmiP 



DmiP 



A vanishing total result at q 2 = is also obtained by Miller, who has studied a family 
of three nonrelativistic constituent quark models. 0] As well, his work suggests that the 
g 2 -dependence is very mild: 

GXf(0) - G^(-0.25 GeV 2 )] < 0.001. (78) 

Capstick and Robson have work in progress that employs a relativized constituent quark 
model. || Using a light-cone meson-baryon fluctuation model, Ma has reported the fol- 
lowing allowed range |§: 

G u ff(0)] Ma = 0.006 -»■ 0.088. (79) 

In light of the uncertainties assigned to HBChPT, there is no essential disagreement 
between the present work and any of these models. Certainly the tendency of nonrelativis- 
tic quark models to prefer G^f(0) = is not obtained from the pion cloud in HBChPT, 
but recall that the effect of higher-energy physics remains unspecified in HBChPT. 
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In conclusion, the effects of isospin violation on the proton's neutral weak magnetic 
form factor have been studied up to next-to-leading order in heavy baryon chiral pertur- 
bation theory. The momentum dependence contains no free parameters, and comes solely 
from the neutron-proton mass difference despite the large number of parameters in the 
Lagrangian. Normalization of the isospin-violating contribution at q 2 = is not specified 
by chiral symmetry, but the pion-cloud effects can be extracted and their contribution is 
roughly 0.02 nuclear magnetons. 
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